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I. INTRODUCTION 



The total energy E12 of two static disjoint objects may be decomposed as, 

E12 = Eo + AEi + AE2 + AE12, (1) 

where Eq is the energy of the vacuum (medium) without objects, AEi and AE2 are (self)-energies required to 
create the objects individuahy in isolation, and AE12 is the change in energy due to their interaction. The interaction 
energy AE12 is finite for disjoint objects if it is mediated by an otherwise free quantum field whose interaction with the 
objects is described by local potentials. It is the only contribution to the total energy that depends on the position 
and orientation of both objects and determines the forces between them. Casimir found that the electromagnetic 
force between two parallel neutral metallic plates does not vanish [3] and that the associated Casimir energy AE12 
may be interpreted as arising from changes in the zero-point energy due to boundary conditions imposed on the 
electromagnetic field by the metallic plates. Zero-point energy contributions to the energies AEi of individual objects 
in general diverge but the change AE12 due to the presence of two disjoint objects is finite. 

Reliable extraction of finite Casimir energies for a long time appeared to be restricted to very special geometries, 
like parallelepipeds [l-Q: spheres d, @, and cylinders [J- A multiple scattering formulation for computing Casimir 
energies of smooth objects was developed by Balian and Duplantier Q but it relied heavily on idealized boundary 
conditions. Kenneth and Klich only recently observed flO] that AE12 may be computed independent of single-body 
contributions to the energy and is always finite for disjoint objects. This two-body interaction energy is compactly 
expressed [l^. Ell terms of the free Green's function. Go, and transition operators, Ti, T2, associated with the 
individual objects, 

AE,2 = \f^ Tr In [l - G0T1G0T2] = 1 ^ g Tr In [l - f,f2] , (2) 

where we have defined partly amputed transition operators Ti — G^Ti, i — 1,2. For potential scattering the interaction 
energy may equivalently be expressed (l2| in terms of the potentials Vi and the corresponding Green's functions Gi 
satisfying GiVi = GoT^. 

In deriving Eq. ([2]) one formally subtracts divergent self-energies and avoids the question of whether these divergences 
have any physical significance. One in particular circumvents the issue raised in [i3,[l3| of how they should be treated 
in the context of gravity, a problem that has so far only been considered for parallel plates [l5l - [T7| . Although it does not 
address such conceptual points, the irreducible contribution of Eq. ([2]) suffices to explain experimental measurements 
of Casimir forces between two disjoint objects. Since the interaction energy for disjoint objects is finite, errors due 
to numerical or other approximations to Eq. ^ can be controlled. This has now been demonstrated by explicit 
calculations for a number of geometries and physical situations p^l - l2l| . 

In this article we examine a recently proposed extension of these ideas to more than two bodies. It was shown in 
[2^ that the irreducible A'^-body part of the total energy remains finite if the N objects have no common intersection. 
We here explicitly evaluate the irreducible three-body part, A_Ei23, of the total energy, 

E123 = Eo + AEi + AE2 + AE3 + AE12 + AE23 + AE31 + AE123, (3) 

in several cases and verify that AE123 remains finite even as irreducible two-body contributions diverge. The formal 
expression for AE123 given in [g^l is considerably more involved than Eq.©, 

A£'i23 = - j 2^ '^^^'^ ^ + Xi2^lT2T3X23 + T1T3T2X32 — TiT2T^ , (4) 

where the X^ 's are solutions to the integral equations, 

[l~t,f,]X,,^l. (5) 

We here obtain and evaluate alternate expressions for irreducible Casimir energies that do not involve a logarithm. 
The article is organized as follows. In Section [ll] we derive Faddeev-like equations for the scattering matrix associated 
with A''-objects and extract their irreducible iV-body parts. The associated A^-body Green's functions are expressed in 
terms of one-body transition matrices describing scattering off each object individually. The procedure is illustrated 
for N = 2, and iV = 3, for which explicit solutions are obtained. The method is generalizable to higher N. The 
general solutions of Section |TT] are used to obtain the Green's functions for two and three semitransparent plates in 
Section [ml The irreducible three-body contribution to the Green's function for three semitransparent plates is found 
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to exactly cancel the two-body interaction of the outer plates when Dirichlet boundary conditions are imposed on the 
central plate. 

In Section HVl we express the irreducible iV-body contribution to the Casimir energy in terms of the A^-body part of 
the transition matrix. This avoids the computation of the logarithm of an integral operator, but requires one to solve 
a set of linear Faddeev-like integral equations for the A^-body transition matrix. In Section|V]we use this formalism to 
obtain irreducible two-and three-body contributions to the Casimir energy of two and three semitransparent plates. 
The three-body contribution again cancels the two-body interaction from the outer plates when Dirichlet boundary 
conditions are imposed on the central plate. 

In Section IVIl we specialize to the case when two of the three potentials are weak. For point-potentials we prove 
that the irreducible two-body Casimir energy is always negative whereas the irreducible three-body contribution 
is positive. The proof immediately generalizes to any form of the weakly interacting potentials. We also derive 
expressions for irreducible two-and three-body contributions to the Casimir energy when the weak potentials have 
translational symmetry and the third potential represents a Dirichlet plate parallel to the symmetry axis. Some of the 
expressions for the irreducible two-body contributions appear to not have been noted in earlier studies. We obtain the 
irreducible three-body contribution to the Casimir energy in this semiweak approximation, and verify independently 
that it is positive and finite. 

In Section [VIII we use these results to investigate the irreducible three-body Casimir energy of a weakly interacting 
wedge placed atop a Dirichlet plate forming a waveguide of triangular cross-section. The potentials forming the 
triangular waveguide overlap and the irreducible two-body contributions to the vacuum energy diverge. However, 
the irreducible three-body Casimir energy is well defined as long as the supports of the three potentials have no 
common overlap. The irreducible three-body Casimir energy is minimal (and vanishes) when the shorter side of the 
wedge is perpendicular to the Dirichlet plate. Inspired by the study in [2^, we investigate the dependence of the 
irreducible three-body Casimir energy on the cross-sectional area and perimeter of the triangular waveguide. To 
emphasize that the finiteness of the irreducible three-body Casimir energies is not only due to the subtraction of 
corner divergences, we, in Section [Villi consider a waveguide with weakly interacting sides of parabolic cross-section 
that touch a Dirichlet plate. The conclusions for weak triangular-wedges generalize to parabolic-wedges with only 
minor changes in interpretation. 

The explicit calculations in this article support the general results of 22] . The irreducible three-body contribution 
to the Casimir energy in the examples considered here is always positive. It furthermore is continuous (and in this 
sense is analytic in the corresponding parameter) when two of the bodies approach each other and intersect. 



II. MANY-BODY GREEN'S FUNCTIONS 



The free Green's function of a massless scalar field in Euclidean space-time satisfies the partial differential equation 

[-V2 + C']Go(x,x')=<5(3)(x-x'), (6) 

where is the Laplacian of flat three-dimensional space. It is related to the corresponding free Green's function 
of Minkowski space-time by a Euclidean (Wick) rotation. The "one-body" Green's function, d, associated with the 
time- independent potential, T^i(x), describing the interaction with the i-th object satisfies 

[-V2 + e + y,(x)] G,(x,x') = 5(3) (x - x'). (7) 

The two-body Green's function solves a similar equation with the potential (Vi + Vj) associated with a pair of 
objects, Gijk denotes the three-body Green's function to the potential {Vi + Vj -I- Vfe), etc. The potentials Vi{x) of 
this model [13] are proportional to (5-functions that simulate the interaction of the scalar field with classical objects. 
Infinitely strong J-function potentials enforce Dirichlet boundary conditions at the surface of the objects. 

One obtains a formal solution to Eq. ([7|) by considering the differential operator as an integral kernel, and manip- 
ulating the kernels as if they were matrices. To emphasize the correspondence between integral kernels and matrices, 
we replace: [— 4- C^] S^^^{-x. — x') G^^, V"i(x) ^(^^(x — x') Vi, and S^^\-k - x') 1. Using ordinary matrix 
algebra one obtains the formal solution to Eq. ([7]) in the form 

Gi — Go ~ GoTiGo, (8) 

where the transition matrix Ti is given by 



T, = V{1 + GoV)-^ = (1 + VGor'V ^V- VGoV + VGoVGoV (9) 
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The second term in Eq. ([8]) is interpreted as due to scattering off the i-th object. It corresponds to the integral 
operator, 



GoT.Go ^ / d^xi J d^x[ Go(x - xi)r,(xi, x;)Go(x; - x') 



(10) 



In the following, symbolic equations often are more compactly writterQ in terms of partly amputated operators. 
Equations for the physical operators are obtained by replacing every partly amputated Greens- function Gi, potential 
Vi^ and scattering matrix Ti, by. 



Gi — ?> G^Gq ^, Vi — > GoVi, 



and 



GoTi. 



(11) 



A. Many-body scattering theory 

The partly amputated iV-body Green's function satisfies the equation 



I + V1+V2 + 



Vn 



G1...N — 1- 



(12) 



The numbers in the subscript of Gi. jv relate to the respective potentials. We may treat the sum of potentials in 
Eq. (jl2p as a single potential and thus proceed as for a single-body. The solution may again be written in the form 



Gi...Ar — 1 — Ti.jv, 

where the A^-body transition matrix Ti. jv satisfies the equation 



1 + {Vi+V2 + --- + Vn) 



Ti., 



.N 



{Vi+V2 + --- + Vn). 



(13) 



(14) 



The solution to Eq. (fT4|) is an infinite series similar to the one in Eq. ([9|), whose terms can be regrouped into 
components T^'' ^, that begin with the i-th potential and end with the j-th potential. For TV potentials we have 
N"^ such components representing transitions from the i-th to the j-th object. This decomposition of the A^-body 
transition matrix is of the form. 



N N 

T1...N = ^^Tl^...N = Sum[f i...Ar], 

i=i ]=i 



(15) 



where the symbol Sum[A] stands for the sum of all elements of the matrix A. The matrix form of the A^-body 
transition operator is. 



Ti 



.N 



1 ^' ' ' ' 1...N 



12 

N -^l...N 



7^21 7=122 712^ 
^l...N ^l...N '■' ^l...N 



\ -'^l...N -'^l...N 



(16) 



where each component is an integral operator. 

Inserting Eq. ([15]), and introducing Kronecker-5 integral operators, Eq. ([T4)) is equivalent to the following set of 
integral equations 



J2 



f^^ — ■ 



(17) 



In matrix notation this set of equations is 

[1 + Vdiag + 01^..Ar] • T1...JV = Vdiag, 



(18) 



^ This is like setting Go = 1- 
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where we have introduced general matrix symbols 





( ^ 


Ax 


Ai ■ 


■ A,\ 




/A, 





• 


■ ° \ 









A2 ■ 


■ A2 







A2 


• 


• 


^l...N — 




A3 


• 


■ A3 


J -^diag — 








A3 ■ 


• 






An 


An ■ 


■ ) 




I 





• 


• An) 



Using these definitions with Eq. ([9]) we write, 

[1 + Vdiag] • Tdiag = Vdiag and [1 + Vdiag] • &I...N = ®1...N> 

and use in Eq. (fT8|) to derive 



.N 



diag- 



(19) 



(20) 
(21) 



The set of linear integral equations of Eq. (|2T|) are often referred to as Faddeev's equations [25|, |26| for nuclear many- 
body scattering, but apparently have been known [l^, [13 in the context of "optical models" for atomic nuclei since 
the 1950's and may have been used in earlier optical studies. The closely related approach in [2^, [s^] is known as 
Martin-Schwinger-Puff many-body theory. 

Eq. (fT3| relates the A^-body Green's function Gi,,,,n to the A^-body transition matrix Ti....n satisfying Eq. (|14|) . 
Faddeev's equations of Eq. (I^TI) reduce the problem of solving Eq. (ITil) for the A^-body transition matrix to that of 
inverting [l-|-0![" ^] by solving a set of N linear integral equations. Remarkably, 0^ ^ depends only on single-body 
transition operators. The norm of 0^ ^ is less than unity (because the norm of single-body transition matrices is) 
and Faddeev's equations can, at least in principle, be solved by (numerical) iteration [26[. 



B. N — 2: Two-body interaction 

Using Eq. (fT2|) the Green's function equation for = 2 has solution given by Eq. (fT3| , where the transition matrix 
is obtained by inverting the Faddeev's equation in Eq. (|2T|) to yield 



ti2 = [1 + 0^2] ' • t 



diag 



Xl2 

X21 



Ti -T1T2 



(22) 



The integral operators in Eq. ((22)) satisfy Eq. Summing the components of T12 we obtain the total transition 
matrix as, 

fi2 = Sum[ti2] = [1 - X12G1] + [1 - X21G2] . (23) 
The total two-body transition matrix T12 can be decomposed into its irreducible one-and two-body parts, 

fi2=ti+t2+ATi2, (24) 

with 



Ti 






f 2 



AT 



12 



X12 
X21 



T1T2T1 ~TiT2 

-f2fi f2fif2 



(25) 



The irreducible two-body transition matrix AT12 includes all contribution with scattering off both potentials. Eq. (j24p 
and the definition in Eq. (1131) imply the following decomposition of the partly amputated two-body Green's function 



G12 = 1 - Ti - r2 - ATi: 



(26) 



where Afi2 — Sum[ATi2] . Summing the four independent two-body transitions of Eq. (|25p . the irreducible two-body 
transition operator in terms of the single-body transition matrices Ti and T2 is. 



Afi2 = Sum[Af 12] = (1 - Xi2)Gi + {I - X2i)G2. 



(27) 
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C. A'^ = 3: Three-body interaction 



One proceeds similarly for three bodies. In this case the formal solution to the Faddeev's equation, of Eq. (|2T|) . is 



- 123 



1[23] 











-'^2 [31] 











-'^3[12] 



Ti —T1G3T2X32 -~TiG2T:iX23 

'T2G3T1X31 T2 -~T2GiT:iXi3 

-T3G2T1X21 —T3G1T2X12 T3 



(28) 



AT12 + AT23 + AT31 



where the G^'s are related to T^'s by Eq. (|5]) and the two-body effective Green's functions, Xij, (i ^ j), solve Eq. ([5]). 
The three-body effective Green's functions, X^y/^-^, {i ^ i ^ k)i satisfy the equation 

^i[]k] [1 ~ TiTjk\ — [1 — TiGjTkXjk — TiGkTjXkj] = 1. (29) 

The total transition matrix in this case is 

fi23 = Sum[f 123] = [1 - Xi[23]Gi] + [1 - X2[31]G2] + [1 - X3[i2]G3] . (30) 

The transition matrix may again be decomposed into irreducible parts 

ti23 = ti + ta + T3 + AT12 + At23 + At3i + AT123, (31) 

where 

-'^i2rLr2T'i + X13T1T3T1 —X12T1T2 —X13T1T3 

—X21T2T1 X21T2T1T2 + X23T2T3T2 —X23T2T3 (32) 

— -^3121371 —X32T3T2 X31T3T1T3 + X32T3T2T3 _ 

is obtained using the N ~ 2 expressions of Eq. (I^S)) . The new irreducible three-body part is, 

(1 — X\2 — X\3 + ^1[23])2^1 \T\X2\ — Xx[23\T\G3X2-^T2 [T1X31 — ^l[23]2^lG2^32]73 
AT123 = [T2X12 — X2[3i]T2G3Xi3]Ti (1 — X23 ^ -^^21 + -'^2(31]) 2^2 [T'2-'^32 ^ -'^2[31]2l2GiX3i]r3 . (33) 
. [23X13 — X3[12]73G2Xl2]ri [73X23 — X3[i2]T3GiX2l]r2 (1 — X31 — X32 + X3[i2])T3 _ 

The decomposition of Eq. (|3ip carries over to the decomposition of the Green's function 

G123 = 1 - fi - fa - f3 - Afi2 - Af'23 - Afgi - Afi23. (34) 

Summing the nine independent three-body transitions in Eq. p3p we find that, 

Ari23 — Sum[ATi23] = — (1 — X12 — X13 +Xi[23])G'i — (1 — X23 — X21 +X2[3i])G2 — (1 — X31 — X32 -|-X3[i2])G3. (35) 

Although not quite as obvious as for two-body scattering, closer inspection reveals that each component of AT123 
indeed involves scattering off all three bodies. A similar procedure can be used to obtain scattering matrices and their 
irreducible parts for more than three bodies. 



III. GREEN'S FUNCTIONS FOR PARALLEL SEMITRANSPARENT (5-PLATES 



We now apply this formalism to derive the Green's functions for parallel semitransparent plates of infinite extent 
described by 5-function potentials 

y,(x) = A,,5(z-aO, (36) 

where specifies the position of the i-th plate on the z-axis, and Ai > is the coupling parameter. In the limit 
Ai — )■ 00 the potential of Eq. p6p simulates a plate with Dirichlet boundary conditions. The translation symmetry in 
the x-y plane can be exploited and Eq. ([7|) written in terms of the dimensionally reduced Green's function, ^^(z, z'), 
defined by 

G.(x,x';C)= 1 ^e^'^^-^^-^^-g^i.z.z'-K), (37) 

where xj^ is the component of x in the x-y plane, is the corresponding Fourier component, and — Cl^ + k^, 
k^ = /c^ -I- fcy. Since the potentials of Eq. (jSH) do not depend on the transverse dimensions, the Green's functions 
G1...N for N parallel plates also correspond to dimensionally reduced gi...N- 



A. A'^ = 1: Green's function for a single semitransparent plate 



When substituted in Eq. ([7]), Eq. ([37]) implies that gi{z, z') solves a one-dimensional ordinary inhomogeneous second 
order differential equation with a (5-function potential that can be solved explicitly, to obtain, 

= i- e-'^l-^'l - e-^l^-^le-^^l--!, (38) 

where Ai = \i/2K, and k was defined after Eq. ([37]) . We also arrive at this solution by starting from Eq. ([8]), which 
for the dimensionally reduced Green's function reads 

9i{z, z') = go{z- z') - ri{z)Un{z'), (39) 

where 50(2, z') is the dimensionally reduced free Green's function, and ri{z) — go{z — a.;). Eq. (jG)) implies that 

g,{z,z')^l-e-^\^-^\ r,(,) ^ l-e-\^-< (40) 

It will be convenient to define, 

f,{z) = 2Kri{z) = e^'^l^^'^'l, and fy = e-'^^"''"^^ = 6"'^"'^ (41) 

where for notational convenience we have defined a^- = ja^ — aj|. The dimensionally reduced transition matrix in 
Eq. (p9| is found by summing the series in Eq. ([9]). Translational invariance in transverse directions and the 5-function 
potential render all integrals trivial and the series can be re-summed to give 

Uiz,z')=2KU6iz-a,)S{z'-a,), U = -^. (42) 

In the Dirichlet limit (A^ — >■ 00) the transition matrix simplifies further to tf = 1. Inserting the dimensionally reduced 
transition matrix of Eq. (j42|) in Eq. (pQ)) reproduces the explicit solution of Eq. ([38]) . 



B. N = 2: Green's function for parallel semitransparent plates 



The previous procedure is readily extended to compute the Green's function of N semitransparent plates located 
at z = ai, i = 1,2, ... ,A^, and described by potentials of the form given by Eq. p6p with associated strengths A^. 
Generalization of Eq. p9p in particular gives the relation 



gi...N{z, z') = go{z - z') - r(z)^ 



.N ■ r(z'), 



(43) 



between the dimensionally reduced Green's function gi,,,N{z^ z') and the corresponding components of the dimension- 
ally reduced transition matrix. The vector r(z) constructed out of the dimensionally reduced free Green's function 
originating or ending at a plate is given by 



r(z)^ = [ri(z),r2(z). 



,rN{z)] = 7^ 
Zk 



-K\z-ai\ -K.\z-a2\ 



— K\z — a]s} I 



(44) 



An advantage of this approach is that the Faddeev integral equations, Eq. (PT|) . collapse to algebraic equations for 
the dimensionally reduced transition matrix ti. jv due to the translational symmetry and the (5-function potentials. 
The transition matrix decouples from the r-vector, which leads to considerable simplification in the evaluation of the 
Green's function. A similar simplification occurs for concentric cylinders and concentric spheres. 

The dimensionally reduced two-body transition matrix can be read out from Eq. (j22p once the corresponding Xij 
has been evaluated. With the single-body transition matrices of Eq. (|42|) all integrals evaluate trivially and the 
solution of Eq. ([5]) for Xij is 



X -X ± 



X, = 1 



i 2 tj G 



(45) 



Using Eq. (P^ . the dimensionally reduced transition matrix for two plates is, 



8 



FIG. 1. Three parallel plates. Plates i and k axe separated by distances aij and ajk from the center plate j. 
Eq. (|46|) inserted in Eq. (|43p gives the Green's function for two semitransparent parallel plates. In the Dirichlet limit, 



ti tf = 1, we have A^- — = (1 — e ^'^"'j), and the transition matrix for two Dirichlet plates simplifies to, 



sinh KQi 



(47) 



From Eq. (j24p we similarly obtain the irreducible two-body part of the dimensionally reduced transition matrix as 



^^t^^' - — - 2 AC 



f ■ JJLf. f . 

^3 An '"I- ^3 



2k 

a" 



tiTijtjVjiti tiVijtj 



(48) 



which in the Dirichlet limit simplifies to 



At 



D 



sinh / 



-1 g-K1ij 



(49) 



The two-plate Green's function has been obtained previously [18j in a more direct manner. We reproduced it using 
the multiple scattering method because this approach readily generalizes to more than two plates. 



C. Af = 3: Green's function for three parallel semitransparent plates 



The three semitransparent plates i, j, and fc, of infinite extent and parallel to the xy-plane are described by potentials 
of the form given in Eq. Without loss of generality we assume that < aj < au (see FIG. [T]). In the previously 

introduced notation this implies that a.y -I- ajk = aik- The vector r{z) in Eq. (j44p now has three components. 

The dimensionally reduced three-body transition matrix is obtained by solving Eq. (j29l) for X^j^] using the Xij of 
Eq. p5)) . For three semitransparent parallel plates one finds that, 



Xiijk]Xjk — 
Using Eq. ^ in Eq 



1 



A 

we obtain 



^ijk — 1 iif'ijijl^ji tj'l^jkik^kj ik^ki^i^ik ~^ '^^i'^ij^j'^jk^k'^ki 



(50) 



2k 



~^j'^ji[k] 



ii'^ij[k] ^3 ti'^ik[j] ^k 
^ki ^j'^jk[i] ^k 

_ ikf^ki[j]^i ^kT^kjlij^j ^k^ij 



(51) 



where 



rij[k] = Tij - riktkTkj- 



(52) 
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It is apparent from the solutions for the N = 2 and N = 3 case that terms contributing to the transition matrix in the 
multiple scattering expansion depend exponentially on the length of the path of propagation. Expanding the inverse 
determinants Ajj and Aijk gives all paths contributing to the scattering operator. This is particularly transparent in 
our essentially 1-dimensional example. 

From Eq. p3p the dimensionally reduced, irreducible, three-body transition matrix is similarly evaluated as, 



At 



ijk 



2k 



:i \Aji Aijk j ' J\ Aji Ajk Aijfc J ^ \Ajk Aijk j 



(53) 



It is interesting to consider the situation when Dirichlet boundary conditions hold on the j-th plate between the 
other two plates. Taking the limit tj tj' — 1, the determinant for three parallel plates is found to factorize into a 
product of two-body determinants, 

A.oofc = A.ooAfcoo = (1 - U e-2-«..)(i _ tk (54) 

where replacing the subscript of a plate by 00 denotes Dirichlet boundary conditions on that plate. In this situation, 
r^k[oo] = (1 - tf)f,k = 0, fjoo[fc] = e'"''^ Afcoo, ^'fcooH = 6"°^*= A,,oo, and Eq. ^ simplifies to 



t 



-t 



1 

Ai^ 



* A.o 
]_ 



-tk; 





tkT^ 



This leads to the observation that 



tfeo 



tf +tfe + At,, 



(55) 



(56) 



Comparing Eq. (j56p with the decomposition of the three-body transition matrix into irreducible one-and two-body 
parts in Eq. (PT|) this implies 

At.oofc + At.fc = 0, (57) 

which confirms the notion that modes in the two half-spaces on either side of a Dirichlet plate are independent and 
that correlations between them must vanish. The irreducible three-body correlations in this limit therefore must 
cancel irreducible two-body correlations between objects on opposite sides of the plate. Taking the Dirichlet limit on 
the central plate in Eq. ((53|) this is verified explicitly. 



At,; 



Aik 



At 



ik 1 



(58) 



where we have used Eq. 

Let us finally consider the case when Dirichlet boundary conditions are imposed on all three plates. The three-body 
determinant again factorizes, Af.f. = Af^A£ = (1 - e"^''°*0(l - e~^''"j'=), and 



t^ 



2k 



1 



2 sinh Ktti 



1 



2 sinh KQi 



2 sinh KQi 



2 sinh KG 

+ 





1 



2 sinh KQjk 
1 



2 sinh KQjk 



2 sinh 



2 sinh 



KQjk _ 



At^ + Atf„ 



(59) 



in the limit of three Dirichlet plates. This implies At^^, + At^ = 0, and is explicitly verified by Eq. (l53l) or Eq. (|58p . 



sinh 



ticiik 







-1 





6^'""''= 



-Atf„ 



(60) 



using Eq. dig]) . 
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IV. MANY-BODY CASIMIR ENERGIES 



Casimir energies are finite parts of the vacuum energy that describe its dependence on configurations of macroscopic 
objects. The interaction of classical objects with quantized fields at low energies can be described by background 
potentials. It is known [iJ: HI, l3ll - [33| that such a semiclassical description for the interaction with a quantized field 
suffers of (local) ultraviolet divergences. A proper treatment of the interaction at high energies requires modeling of 
the quantum fiuctuations associated with the objects. One fortunately sometimes is able to isolate parts of the vacuum 
energy that depend only on global changes of the system and can be reliably computed in semiclassical approximation. 
In the following we systematically determine irreducible parts of the vacuum energy for a given number of classical 
objects. These irreducible iV-body Casimir energies diverge only if all N potentials describing the classical objects 
have a region of common support [2^ . 

Let Eo be the (infinite) vacuum energy associated with zero-point fiuctuations of a massless scalar field in the 
absence of all background potentials, Vi{'x). The change in vacuum energy in the presence of N objects associated 
with the background potential, V = Vi, can be derived using field theory techniques, for example in [18, 24], to be 

Ei...N-Eo = -l r ^2C^Tr(Gi...Ar-Go) = -i T ^ Tr In Gi...jv - H ^ Trln(l ~ fi...^). (61) 



2 J-oo 27r 2 ./-^ 27r 2 2n 



These expressions have recently been dubbed the Trace-G-formula and Trace-Log-G-formula respectively. For fre- 
quency independent potentials, the relation between them is established by differentiating Eq. ([7]), 



_d 



^G = GG, (62) 



and ignoring a boundary term. 

To proceed further we generalize Eqs. ((25)) and ([M)) and decompose a Green's function involving N potentials into 
irreducible A^-body parts, 

i i<j 

Using Eq. (I62[) . the irreducible one-two-and three-body parts of the Green's functions can be written in the form 

AG, = G, - Go = --4 In ^, (64a) 
dC Go 



AG,, = G,, - AG, - AG, = In (64b) 

GijkGiGjGk 



AGyfe = G,,fe - AG,, - AG,k - AG^, - AG, - AG, - AG,. = -—^ In ^'-'" ' ^ , (64c) 

which is a (cascading) recursive definition that can be extended to higher N . 

Eqs. (j6ip and (j63p imply a corresponding decomposition of the vacuum energy into irreducible iV-body contributions. 



i 

where 



Ai?i...Ar = -i r ^2C2TrAGi...^. (66) 
2 ./-ex, 27r 



As shown in ^25], and as will be explicitly verified in examples below, the irreducible A^-body contribution to the 
vacuum energy diverges only if all N potentials have a common support. One-body vacuum energies thus are generi- 
cally divergent, whereas two-body Casimir energies diverge only if the two bodies intersect (and thus could be viewed 
as one). More interestingly though, three-body Casimir energies diverge only when all three objects have a common 
intersection-the three bodies need not be mutually disjoint and could, for instance, be arranged to form a triangle. 

Eq. ([T^ relates the irreducible iV-body contribution of the Green's functions to the irreducible A^-body transition 
matrix, 

Tr AG1...JV = -Tr Ari...jvGoGo = Tr Ari...w^Go. (67) 
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The support of delta-function potentials Vi is restricted to the surface Si of the i-th object and components of the 
transition-matrix at most have support on the union of two such surfaces. It is therefore convenient to formally define 
a vector R(x), and a matrix R, with components 



Go(x 



R,(x) :=Go(x-y 
Using these definitions in Eq. (|62l) we have 

-^R = J d32;R(x).R(x) 



yes- 



(68) 



(69) 



Writing the irreducible A^-body transition operator in Eq. (j67p in matrix notation and using Eq. (j69p , we express the 
irreducible A^-body contribution to the vacuum energy of Eq. (j66p in the form 



1 



, zvr 



dC 



AT.... . -^R 



(70) 



The trace in the last expression is over matrix indices and includes integrals over the lower dimensional surfaces of the 
associated objects. Note also that Eq. ([70]) involves the irreducible iV-body transition matrix, Ti...Ar, not its partly 
amputated cousin Ti...jv. 



CASIMIR ENERGIES FOR PARALLEL SEMITRANSPARENT J-PLATES 



We illustrate this formalism by evaluating the irreducible (scalar) Casimir energy for semitransparent parallel plates 
described by potentials of the form given in Eq. (j36]). Exploiting translational invariance parallel to the plates in 
Eq. (j70p , the irreducible TV-body Casimir energy per unit area is described by dimensionally reduced quantities 



N 



1 

6^ 



K dn Tr 



At,... • rf^r 



(71) 



where and Ly are the (infinite) lengths of the plates in x and y direction, respectively. The integrals on fcj,, 
and ky, are performed using polar varibales, which effectively amounts to replacing — >■ 2k^/3, where k was defined 
after Eq. (1571) . The dimensionally reduced transition matrices, Ati...., for = 2 and = 3 are, respectively, given 
by Eqs. (pS)) and (1551) . The derivative of the dimensionally reduced free Green's funct 



ion in this case is the matrix 



1 

(1 -f- Ka2i) e 



(1 + KflATi) e" 



(l + Kai2)e-'"''i= 

-na2i 



(1 + KaN2) e 



{I + KaiN) e-""'" 
(1 -1- Kflajv) e-«"2" 



(72) 



where aij is the distance between the z-th and j-th parallel plate defined previously. 



A. A'^ = 1,2: Irreducible one— and two-body Casimir energy for semitransparent plates 



The irreducible one-body vacuum energy per unit area associated with the i-th plate diverges. Eq. ([7T|) gives it as 
the integral 



LxLy 



1 



127r2 



(73) 



with ti defined in Eq. (|42p . The one-body vacuum energies are ultra-violet divergent at any non- vanishing coupling, 
but do not depend on the relative position of the plates and therefore do not contribute to forces between them. 

The irreducible two-body Casimir energy per unit area associated with plates i and j is obtained by inserting 
Eqs. (gS]) and dZSl) (for A^ = 2) in Eq. 



AE^ 

LxLy 



127r2 





1 

1 


/ K^dK 


Jo 


A,, 



(74) 
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where the two-body determinant is given by Eq. (|45l) . Eq. (|74l) for the Casimir interaction energy of two semitrans- 
parent plates was obtained previously in 2J]. In the Dirichlet limit, — >■ 1, Eq. ([7^ simplifies to the well known 
Casimir energy for a massless scalar field satisfying Dirichlet boundary conditions on a pair of parallel plates, 



LxLy 



1 



127r2 



' K^dK ^'^"'^ 



1 



1440 a3- ■ 



(75) 



Eqs. (f74l) and (|75|) are finite and negative for two disjoint plates. 



B. N = 3: Three-body Casimir energy for three parallel plates 



The irreducible three-body Casimir energy of three plates is similarly obtained by inserting Eqs. (j53p and (j72p (for 
iV = 3) in Eq. ([71]), 



AE, 



ijk 



LxLy 



U<1 



+tk < 1 



1 


1 


Ajk 1 




A^ ~^ 


Aijk J 


1 


1 


^ A,k 




'ATk 


^ijk 


1 


1 




Aki 







2(1 + najk) \ - 1 [> <{ 1 - rjk[i]e 



l-f. 



k\j\e 



+ 2(1 -I- nai 



1 
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rij[k]t 



Ajk 1 


A„fc J 




\ 


Aijk , 




1 


> 


Aijk , 





(76) 



When Dirichlet boundary conditions are imposed on the central j-th plate, the relation between irreducible two- 
and three-body transition matrices noted in Eq. ()57p implies a corresponding relation between two-and three-body 
Casimir energies. 



AE.^k + AE,k = 0. 



(77) 



This is explicitly verified by using the factorization of the three-body determinant in Eq. (|54p and the Dirichlet limits 
for fy[fc] given after Eq. (|54l) in Eq. ([76]) . and identifying the irreducible two-body energy of Eq. ([74|) in the result. 

In the Dirichlet limit for all three plates the irreducible three-body Casimir energy cancels the well-known two-body 
interaction between the outer Dirichlet plates. 



AE^ 

LxLv 



1 



1440 



LxLy 



(78) 



where aik is the distance between the outer plates. 

This cancellation is to be expected on physical grounds and serves to check the calculation. For semitransparent 
plates the cancellation is not complete and the irreducible three-body contribution to the total Casimir energy can 
be significant for parallel plates. Note that the sign of the irreducible iV-body contribution to the scalar Casimir 
energy alternates. Although not apparent from the expression of Eq. (j76p . this irreducible three-body contribution to 
the Casimir energy is positive for any positive couplings Ai, A2, A3 and any relative position of the three plates. For 
parallel semitransparent plates we thus verify the more general result obtained in [22]. Also, as discussed in [2^ and 
noted previously, the three-plate Casimir energy diverges only if all three plates coincide. 

In the following we will see that these generic results for the sign and analyticity of the three-body scalar Casimir 
energy hold in the limit where two of the three potentials are weak and need only be accounted for to leading order. 

VI. THREE-BODY SCALAR CASIMIR INTERACTION FOR SEMIWEAK COUPLING 



We now consider irreducible vacuum energies for three bodies when two of the three potentials, Vi and V2, are weak 
and need only be taken to leading order. No restriction is imposed on the potential V3 describing the third body. To 
the leading order we thus approximate Ti ^ Vi and T2 ~ V2 in Eq. ^ . The three-body transition matrix of Eq. 
in this semiweak approximation simplifies to 



-■- 123 



1 
1 



^3^12] 



Vi -1^1(1-^3)^2 -Viil-V2 + nV2)T3 

'V2{l~f3)Vi V2 ~V2{l~Vi+f3Vi)f3 

-^3(1- 1/2)14 -T3(l-Vi)V2 T3 



(79) 
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Here -^^2] satisfies Eq. ([29)) . which to leading semiweak approximation is solved by 

XY[12\ - 1 + ^3(1 - Vi)V2 + fail - V2)Vi + f3Vif3V2 + f-iV^nVi. 



(80) 



The transition matrix in semiweak approximation of Eq. (j79p may again be decomposed into its irreducible one-two- 
and three-body parts, leading to the semiweak version of Eq. ([3T|) . From Eq. p5|) the irreducible two-body transition 
matrices in semiweak approximation are, 



-V1V2 
-V2V1 







(81) 



with z = 1, 2. Similarly approximating Eq. p3p . the three-body transition matrix in semiweak approximation becomes, 

(82) 



Vinv2 

V2f:iVi V^2G3Fif3 

T3V'2G3^l f3FiG3V'2 -73^103^2^3 - 1-3^2(53 ^lf3 



where G3 = 1 - 73. 

Casimir energies in the semiweak approximation are obtained using Eqs. (j66p and (j67p . Inserting Eq. (j8ip in 
Eq. ((S7)) we have to this approximation. 



TrAG]^ 


= Tr 


AT]^ 


d 

dC 


R 


d 

' dC 


TrAGS' 


= Tr 


AT,^ 


d 

dC 


R 


d 

" dC 



rTr 



G0V1G0V2 
GoViGoT3 



(z = l,2). 



The corresponding irreducible three-body contribution using Eq. ([5^ in Eq. (|S7)l is 



TrAG]^3 = Tr 



ATr23 • 



= -dc^' 



GQViGQT'iGQV2 + GoV2GoT'3GoVi — GoT'3GoViGoT'3GoV2 



(83a) 
(83b) 

(84) 



Inserting Eq. (|83p in Eq. (|66l) . and integrating by parts, the irreducible two-body Casimir energies in semiweak 
approximation are 



^EY2^~\ P ^Tr[GoViGoy2]. 



dC 
2tt 



Tr[GoV^GoT3 



(« = 1,2), 



(85a) 
(85b) 



verifying results reported in |13 |. The corresponding irreducible three-body contribution to the Casimir energy in 
semiweak approximation using Eq. (j84l) in Eq. (j66p is 



A-E'];^3 



1 



f°° 

J ^ Tr[GoViGor3GoV2 -I- G0V2G0T3G0V1 — G0T3G0V1G0T3G0V2] ■ 



2 J -00 2^ 

In the following we evaluate Eqs. (|55|) and for some special cases. 



(86) 



A. Weak point potentials 

Weak point potentials of the form, 

Vi{-x.) = Ai^(x- Xi), 
for i — 1,2, allow one to explicitly perform the integrals in Eqs. (1851) and 

Ai A2 



In this case we have 



AE- 



w 



■/°° §[Go(xi-X2)]'<0, 



(87) 
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and, using Eq. © in Eq. ((85b|) . 

^E^^~Y /3^^{Go(0)-G3(x„x,)}<0, (z = 1,2), (89) 

because the integrand in braces is positive for positive V3. The irreducible two-body contributions to the vacuum 
energy thus are negative for weak point potentials. We similarly obtain that the irreducible three-body correction to 
the vacuum energy, 

A<3 = ^/^|^{[Go(xi-x,)]^-[G3(xi,X2)]^}>0, (90) 

in this case is positive for any (positive) potential V3. Note that the irreducible three-body Casimir energy in semiweak 
approximation diverges only if Xi = X2 is in the support of V3. 

The pattern in the sign of the irreducible A'^-body contribution is consistent with the findings of [22| . Furthermore, 
since any positive potential is a (positive) superposition of point potentials, this pattern of the signs of irreducible 
A'^-body contributions extend to any shape of the objects in semiweak approximation. This is explicitly verified by 
the following examples. 



B. Weak potentials with translational symmetry parallel to a Dirichlet plate 



We consider a Dirichlet plate and weak potentials that do not depend on the Cartesian coordinate x, 

V^i—Vi{y,z), for i — 1,2; and V3 — X3 S{z — a^), with A3 — ^ cx). (91) 

To evaluate Eqs. (|85l) and ([55]) for such potentials we require the operator GoT-^Gq for a Dirichlet plate. In order to 
exploit the translational symmetry in a;-direction we write the solution to Eq. (|6]) for the free Green's function in the 
form 



Go(|xi-x2|;C)- TTT^'^'''-''''"'''^^^ = / -^e^'^==^---)^^^-^ = — (92) 

{2n)^ 2k 2tt 2tt 47r|xi-X2| 



where di2 = \/ (j/i — J/2)^ + (zi — z-i)^ is the projected distance in the x\ = X2 plane, and — k"^ + Ko{x) is 
the modified Bessel function of zero order. Note that k defined after Eq. ([57]) satisfies ^ k? + ky. Using the first 
equality of Eq. (1^^ and the dimcnsionally reduced transition matrix of a Dirichlet plate given in Eq. (H^ one can 
show that 

-AGf(xi,X2;C) = [Gor3^Go](xi,X2;C) = G„(|xi-X2|;C)- f ^ ,^kA^^~^.)^^o{^^ (93^ 

where X2 — {x2,y2, + 203), and di2 is the length of the shortest path between xi and X2 in the {xi — a;2)-plane 
that reflects off the Dirichlet plate. For a Dirichlet plate at z = 03, this distance is given by di2 = (j/i — J/2)^ + (ki — 
flsl + 1^2 — flal)^- A geometrical interpretation of di2 and di2 is shown in FIG. [H Substituting Eq. in Eq. (|5]) 
leads to 

Gf (xi , X2 ; C) = Go ( |xi - X2 1 ; C) - Go ( I xi - X2 1 ; C) , (94) 

which is anti-symmetric under reflection about the Dirichlet plate. Note that if xi and X2 are on opposite sides of 
the plate, di2 = di2, and Gf vanishes. 



Substituting Eq. ((92| for the free Green's functions, and Eq. ((93|) for the irreducible Green's function of a Dirichlet 
plate in Eqs. ([55]) and using the identity 



f 

Jo 



• _ _ _ 1 

KdKKo{Kx) = —, (95) 
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FIG. 2. The distances di2 and di2. The effective distance di2 is the shortest distance between the two points for a path that 
reflects off the Dirichlet plate at z = aa. It also is the shortest distance between {yi, zi) and a mirror image of the point (1/2, Z2) 
with respect to the z — line. 



the irreducible two-body Casimir energies per unit length in semiweak approximation for potentials with translational 
symmetry are 



Lr 327r2 



{i^l,2). (96b) 



The Casimir energy in Eq. (|96ap for two weakly interacting objects with translational symmetry was previously 
obtained in [3^ . The Casimir energy for a Dirichlet plate weakly interacting with an object possessing translational 
symmetry was obtained in [l2;j, but was given as a series involving modified Bessel functions. The expression in 
[l2| generally is much harder to evaluate than Eq. (j96b|) . The simplification in Eq. (j96b|) was achieved by using the 
effective Green's function for a Dirichlet plate in Eq. (|93|) . For many potentials, the evaluation of the Casimir energy 
by Eq. (I96bp is immediate. We can for example calculate the two-body Casimir energy for a cylinder of radius a, 
described by the weak potential V^yi — Xd{r — a), interacting with a Dirichlet plate positioned sd z = R > a. From 
Eq. ((96bl) one readily finds, 



327r2 7o 7o |rsin6i-i?|2 16tt 



i?2 



(97) 



which reproduces the expression in |12| . A similarly simplified evaluation is expected for an arbitrary surface with 
translational symmetry weakly interacting with a Dirichlet plate parallel to the symmetry axis. 

The irreducible three-body Casimir energies for translationally invariant weak potentials and a Dirichlet plate can 
be similarly evaluated using Eq. (1551) . The first two terms in Eq. involve the product of the free Green's function. 
Go, with the irreducible Green's function for a Dirichlet plate given in Eq. (|93|) . The last term requires the product 
of two irreducible one-body Green's functions. A useful identity for the product of two modified Bessel functions of 
zeroth order is 

' 1 / \ 1 

RdRKo{Rx)Ko{Ry)^- jln - (98) 



Inserting Eqs. (|92|) and (|93|) in Eq. (|86)) to write the Green's functions in terms of modified Bessel functions, and then 
using Eq. (jM]) . we obtain 

where the distances di2 and di2 were introduced earlier and are shown in FIG. [H The function 



(100) 
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{0,h) 




FIG. 3. Weakly interacting triangular-wedge on a Dirichlet plate. The objects are of infinite extent in the x-direction. The 
weakly interacting sides of the wedge (in red) have finite length. 

is bounded by 1 < Q(x) < 1 — 21na:: in the relevant domain < a; = < 1. This implies that the three-body 

"12 

Casimir energy of Eq. (j99p is always positive and bounded by 

di2 is the distance between a point on the first object and another point on the reflected image of the second object 
(see FIG. [2]). It vanishes only at points where the two weak objects and the Dirichlet plate are concurrent. The 
irreducible three-body Casimir energy in the semiweak approximation of Eq. (|99p thus is finite if the three objects 
have no point in common. This contribution in particular does not diverge as the objects approach the plate or each 
other, corroborating the findings in (22]. Note that the lower bound in Eq. (|10ip is the two-body Casimir energy 

between weak potentials of Eq. (|96ap . but with the reflected object (di2 — > ^12) and of opposite sign. The irreducible 

^2 

three-body Casimir energy approaches the lower bound for ~ 1 and thus partially cancels the irreducible two-body 

"12 

energy if one or both objects approach the Dirichlet plate. In fact, if the two weakly interacting objects are entirely 
on opposite sides of the Dirichlet plate, the lower bound is achieved because di2 — di2 and the three-body Casimir 
energy cancels the two-body interaction energy between them precisely. 

The following examples demonstrate the finiteness, sign, and analyticity, of three-body contributions to Casimir 
energies for cases in which irreducible one-and two-body contributions to the vacuum energy diverge. 



VII. TRIANGULAR- WEDGE ON A DIRICHLET PLATE 



We first consider a triangular- wedge with two sides described by weak potentials atop a Dirichlet plate at z = 0, 
forming a waveguide of triangular cross-section: 

Vi{y, z) ~ XiS{—z + niaiy ~ a)) 9i, with 9i = 9{y — mm[0, a]) 9{max.[0, a] — y), (102a) 
V2{y,z) = X26{-z + mp{y~b))92, with 92 = 9{y - mm[0,b]) 9{m&x[0,b] - y), (102b) 
Vsiz) = A3^(z), with A3 -> 00. (102c) 

The sides of the wedge have slopes rua = — cot a and = — cot /3, and lengths y/hF+a^ and y/h^ + lP, respectively. 
The constraint iriaa — mpb — —h forces the sides to intersect at {y = 0, z = h), where h is the height of the triangle. 
The base of the triangle formed then measures \b — a\. Note that the Dirichlet plate at z = is of infinite extent. 
This triangular- wedge on a Dirichlet plate is depicted in FIG. [31 Suitable parameters for describing the triangular 
waveguide are {h,a,f5), or (/i, a — a/h,b — b/h). Without loss of generality we measure all lengths in multiples of 
the height h. The triangle then has height h — 1 and the parameter space for the triangle is —00 < a,b < 00, or, 
equivalently, —tt/2 < a,f3 < tt/2. 



Observe that all irreducible two-body Casimir energies in Eq. (j96|) diverge due to ultra-violet contributions from the 
corners of the triangle where pairs of potentials overlap. More precisely, the integrand AE^ diverges when di2 ~ 
near the vertex of the wedge. The integrand of diverges when Zj ~ near the corner with the Dirichlet plate. 

The irreducible three-body Casimir energy, AEj^^, in Eq. (j99p on the other hand is finite because di2 never vanishes 
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FIG. 4. Casimir Landscape: f (a,/3) as a function of the opening angles a and /? for a weakly interacting triangular- wedge 
on a Dirichlet plate. The valley connecting the a > P region with the a < /3 region is an artifact caused by limited numerical 
accuracy. The valley should be replaced by a very thin and infinitely high wall describing the sharp change in energy when all 
surfaces overlap. On the right, the shapes of the triangles are matched to the respective regions of the a-j3 plane. 



in the integration region. Substituting the potentials of Eq. (|102|) for the semiweak triangular waveguide in Eq. (|99|) 
and evaluating the z-integrals gives, 

JO JO "12 \"12/ 

where we have rescaled the integration variables, yi — \a\ui and y2 = \b\u2 by the respective lengths. All distances 
have been expressed in units of h: di2 — hui2 and di2 — hui2, with 

UI2 = {aui - hu2f + [|1 - liil + |1 - ^t2|]^ (104a) 

u\2 = {aui - hu2f + (wi - ^2)^. (104b) 

With the function Q{x) defined in Eq. (|100l) the three-body interaction energy of Eq. (|103l) is finite and can be 
evaluated numerically. In FIG. |4]we plot £(a,/3) as a function of the angles a and /3. The three-body interaction 
energy is always positive and vanishes (and is minimized) only for a = 0, or /? = 0. It is minimal when the shorter 
side of the wedge is perpendicular to the Dirichlet plate. Wedges with angles /3<0<aora<0</3 are energetically 
preferred over wedges with angles a, /3 > or a, /3 < 0. The three-body Casimir-energy diverges only when all three 
sides of the triangle have a point in common, i.e. when a = /3, or a = — /3 = ±7r/2. 



f(a,/3) 



A1A2 

327r3 



Abalo, Milton, and Kaplan, recently |23[ investigated the dependence of the Casimir energy on the area and perime- 
ter of triangular waveguides on which Dirichlet boundary conditions were imposed. Although only interior modes were 
taken into account and divergences associated with corners and single-body vacuum energies were removed ad hoc, 
they found that the dimensionless Casimir energy of their triangular wave guides closely follow a universal function 
of the dimensionless ratio (P'^/A) of the perimeter P and area A of the cross-section. This would imply that the 
Casimir energy of triangular wave guides depends on just one, rather than two, dimensionless parameters. Although 
we cannot expect a similar dependence, the universal behavior observed in i23j prompted us to also investigate the 
dependence of the semiweak three-body Casimir energy on the dimensionless perimeter p = [P/ h) and dimensionless 
area s = {A/h?) of the triangular waveguide. It is also of interest to inquire for what configuration the energy of a 
triangular waveguide is minimized if the cross-sectional area is kept fixed. The dimensionless area s and perimeter p 
of the triangular wedge are given by, 

\\h-dl (105a) 
\h-h\ + ^Jl + ^? + \Jl + y^. (105b) 

The parameter space of a triangular-wedge in this case is s > 0, and p > 2s + 2\/l + > Max(2, 4s). See FIG. H) 
The inequality, p > 4s, is a consequence of the triangle inequality. 



A 

P 
'h 
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I Reflections 

FIG. 5. £{a,b) as a function of a for fixed area, A = h^. The irreducible three-body Casimir energy is minimal when the 
shorter side of the wedge is perpendicular to the Dirichlet plate (5 = or 6 = 0). The maximum in the intermediate region 
corresponds to the unstable equilibrium of an isosceles triangle. The dashed curves are the approximation £{a, b) ~ \ab\ obtained 
by replacing the integrals in Eq. (|103p with unity. The dotted curves are reflections about the £ = line. 



In FIG. Owe plot the energy as a function of a for fixed area: A — h^, oi \b 



2, or I tan /3 - tana| = 2. The 



three-body Casimir energy for a waveguide of given cross-sectional area is minimal when the shorter side of the wedge 
is perpendicular to the Dirichlet plate (a = or /3 = [a = tan~-'^(— 2)]). In the intermediate region the energy is 
extremal for an isosceles triangle (—5 = 6=1) with £{—1, 1) = 0.893112 . . .. The dashed curve in FIG. [5] represents 
the approximation £{d,b) ~ \db\ obtained by setting the dimensionless integral in Eq. (|103|) to 1. Remarkably, this 
extremely simple expression for the irreducible three-body energy is accurate to better than ten percent everywhere. 
We also show reflections of the curves to illustrate that the discontinuities in the slope are entirely due to the absolute 
value in the pre-factor \ab\ and the integral itself is analytic. 

We rewrite the irreducible three-body Casimir energy as a function of the cross-sectional area and perimeter by 
inverting Eqs. pOSp to obtain 




if b > d, 

if 6 < 5, 

if 6 > 5, 

if 6 < a, 



where 



M = 



1 {p-2S) 



p{p - 4s)|p(p - 4s) - 4| 



2p [p - As) 

Substituting Eqs. (|106p in Eq. (I103p . the three-body Casimir energy as of function of perimeter and area is 



£{s,p) 



^ duidu2 



Q 



'0 Jo "12 

where the rescaled distances in terms of area and perimeter are given by 

l2 



'12 



"12 
"12 



[/2(mi - U2) + s{ui + U2)] + [|1 - Ul| + |1 - U2I 
[/2(mi — U2) + s{ui + 112)]^ + (ui - 1*2)^. 



(106a) 
(lG6b) 

(107) 

(108) 



(109a) 
(109b) 



In FIG. [S] the irreducible three-body contribution to the vacuum energy of a semiweak wedge is plotted as a function 
of dimensionless area and perimeter of the cross-section. The energy now is minimal along the curve 



1 + 2S + Vl + 4|2 



2 + 2.5 + 0(25)2 
l + 45 + 0(^) 



2S < 1, 
2S > 1, 



(110) 
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which corresponds to right-angled triangles. The energy diverges along the line p > 2, s = 0, which corresponds to the 
two sides of the wedge coinciding (a = /3). In FIG. |6] the curve p = 1 + 2s + Vl + 45^ for s > corresponds to right 
triangles of minimal energy, and the boundary of the parameter space at p = 2s + 2y/l + s^ for s > is associated 
with isosceles triangles. 

We do not observe that £{s,p) is a function of p^/s only. The rather good approximation obtained by ignoring the 
dependence on the integral in Eq. (|108p . suggests that the energy approximately is given by 

r2 ~2| mp-2sr-pHp~A~sr\ 

ns,p) - IM - ^ I = ^^^^^ , (111) 

a somewhat involved function of the perimeter and area. 



VIII. PARABOLIC- WEDGE ON A DIRICHLET PLATE 



To explicitly verify that not just corner divergences have been subtracted in the irreducible three-body contribution 
to the vacuum energy [22;], we also consider a weakly interacting parabolic- wedge on a Dirichlet plate. It is described 
by the potentials 

Vi{y,z) = \iS{-z + a{y - a)'^)Oi, with 6i = d{y - mm[0,a]) 9{m'Ax[0, a] - y), (112a) 
V2{y,z) = \2Si-z + p{y-b)^)e2, with 02 = e{y ~ mm[0,b]) 9imax[0,b] ~ y), (112b) 
V3{z) A3(5(z), with A3 00. (112c) 

The parameters a and /3 here give the foci of the parabolas and have dimensions of inverse length. The constraint 
aa^ = ph"^ = h implies that the two parabolas intersect at {y = O^z = h). See FIG. [71 As in the case of the wave 
guide with triangular cross-section, the base has length |6 — a| and the height of the wedge above the Dirichlet plate 
is h. The parameter regions are: —00 < a,b < 00, or, equivalent ly, < a, l3 < od. We measure lengths in multiples 
of h and use parameters {h, a = a/h, b = b/h) to describe it. 



We proceed exactly as for the triangular-wedge and find that the three-body Casimir energy of a parabolic-wedge 
is also given by Eq. (jl03p . except that the distances now are given by 

ul^ = {aui - hu2f + [(1 - uif + (1 - U2fY, (113a) 
UI2 = (5wi - hu2f + [(1 - uif - (1 - U2ff. (113b) 
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(-a,0) (6,0) 
FIG. 7. Weakly interacting parabolic-wedge on a Dirichlet plate. 




— £{a,b) 
-- \ab\ 

— Reflections 



FIG. 8. Irreducible three-body Casimir energy for a parabolic waveguide of given cross-sectional area. See FIG.[5]for description. 



The three-body Casimir energy of a parabolic- wedge on a Dirichlet plate also is minimized when either a — 0, a = oo, 
or 6 = 0,/3 = cx). Due to the constraint, aa"^ = /36^ = h, the shorter side of the parabolic wedge in this case degenerates 
to a straight line perpendicular to the Dirichlet plate. Most of the analysis of the waveguide with two sides of parabolic 
cross-section is the same as for a triangular one with only minor changes in interpretation. We note that the rescaled 
area and perimeter of the parabohc wedge are 



P^\b-a\ + - 



1 + - 

a* 



&2 , 2 / 4 



(114a) 
(114b) 



The three-body Casimir energy of a semiweak parabolic-wedge is shown in FIG. |S1 The approximation of replacing 
the integral by unity is not very accurate in this case, but the overall dependence of the irreducible three-body energy 
of a parabolic waveguide on the parameters a and b is qualitatively similar to that of a triangular one. 



IX. DISCUSSION 



In Casimir studies one generally is interested in dependence of the vacuum energy of massless quantum fields on the 
presence of objects whose interaction with the quantum fields is treated semiclassically, with quantum fluctuations 
of the fields describing the objects themselves being disregarded. This leads to an effective action with ultraviolet 
divergent contributions associated with geometrical properties of the objects reflected by the coefficients [H, [sl] in 
the asymptotic expansion of the heat kernel [sil [stI l38j|. The corresponding ultra-violet divergences in the vacuum 
energy are proportional to the spatial volume, surface areas, curvatures, as well as number and type of corners or 
intersections of the objects. They depend only on local geometric properties of the system. 
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Fortunately, there also are non-local contributions to the vacuum energy that do not depend on the high energy 
description of the model and can be reliably obtained in semiclassical approximation. The best known of these is 
the force between disjoint classical objects due to vacuum fluctuations, first obtained for parallel metallic plates by 
Casimir T|. It has since been shown that this force is always finite [l3| and that the associated finite contribution 
to the vacuum energy may be computed for arbitrary obj ects in terms of the single-body scattering matrices with 
Eq. ([2]). The investigation of generalized pistons in [3^, |40] suggested that one may isolate finite parts of the vacuum 
en ergy that describe the forces between objects even if these are not mutually disjoint. These ideas were formalized 
in |22| where irreducible N-body contributions to the vacuum energy were defined recursively and shown to be finite 
unless the N-bodies have a common intersection. For a scalar field whose interaction with N-objects are semiclassically 
described by positive local potentials, the irreducible contribution to the vacuum energy furthermore was found to be 
positive for an odd, and negative for an even, number of objects. 

We have here put these general considerations on a much more practical and concrete footing and developed a 
formalism to extract and compute irreducible iV-body contributions from the single-body transition matrices. Starting 
from Faddeev's equations in Eq. ([2T|) . the irreducible parts of the TV-body scattering matrix were extracted recursively. 
We used this formalism to compute several examples of irreducible two-and three-body Casimir energies. All our two- 
body results have been obtained previously, but we were able to reproduce some of them in a much simpler and direct 
manner. Our three-body results for irreducible Casimir energies are new. The irreducible three-body contributions 
to the Casimir energy of parallel semitransparent plates was obtained analytically and indeed remains finite when 
two of the three plates overlap. We showed explicitly how the irreducible three-body contribution precisely cancels 
the irreducible two-body Casimir energy of the outer plates when Dirichlet boundary conditions are imposed on the 
central plate-providing a raison-d'etre for both, the existence, and sign, of the three-body contribution to the force. 
For semitransparent plates the cancellation is not complete but can be sizable. 

In Section [VII we analyze the irreducible three-body interaction in semiweak approximation. In this approximation 
we are able to compute the irreducible three-body Casimir energy for objects that are not mutually disjoint and 
whose irreducible two-body contributions diverge. The irreducible three-body contributions to the vacuum energy of 
a waveguide constructed by placing a weakly interacting triangular-or parabolic-wedge on top of a Dirichlet plate 
was found to be finite and computable without intermediate regularization. Our examples demonstrate that not only 
corner divergences, but also divergences related to curvature are subtracted by this procedure. We also explicitly 
verified that the irreducible three-body contribution to the vacuum energy of a massless scalar field is positive. 

To develop a better understanding in a non-perturbative setting, we are currently investigating the irreducible 
three-body vacuum energy of a triangular waveguide formed by imposing Dirichlet boundary conditions on three 
intersecting infinite planes (the geometry is similar to that of FIG. [31 but with sides of infinite extent). In the limit of 
an extremely flat triangular cross-section, we intend to compare the numerical results with analytic calculations. We 
further wish to extend these methods to the physically relevant electromagnetic case. Although irreducible three-body 
contributions to the vacuum energy are expected to remain finite, we so far have no rigorous statements about their 
sign for vector fields. Interestingly, it is at least conceptually feasible to directly measure irreducible electromagnetic 
three-body contributions to the vacuum energy by balancing off irreducible two-body parts. We are investigating 
whether this is experimentally feasible. 
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